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PURPOSE:  The  Coastal  Engineering  Technical  Note  (CHETN)  herein  describes  the  EbbJet 
Calculator,  a  utility  of  the  Analytical  Toolbox  of  the  Diagnostic  Modeling  System  (DMS)  (Kraus 
2000).  The  EbbJet  Calculator  estimates  depth-averaged  velocity  and  sediment  transport  rate  at 
an  inlet  from  established  analytical  and  empirical  equations. 

BACKGROUND:  As  described  in  CHETN  IV-28  (Kraus  2000),  “Introduction  to  the 
Diagnostic  Modeling  System  (DMS),”  the  DMS  is  an  experienced-based  methodology  for 
addressing  sediment-shoaling  problems.  It  comprises  a  suite  of  diagnostic  tools  intended  for 
planning,  design,  and  maintenance  of  navigation  projects.  The  DMS  methodology  provides 
engineering  guides,  procedures,  and  analytical  tools  founded  upon  the  principles  of  fluid 
dynamics  and  sediment  transport.  As  a  diagnostic  procedure,  the  DMS  is  not  intended  to  replace 
the  level  of  detailed  information  required  for  final  design  of  newly  authorized  channels  or  for  in- 
depth  study  of  an  existing  navigation  project.  Rather,  it  is  intended  to  function  as  a  screening 
tool  to  quickly  assess  existing  shoaling  problems,  to  categorize  them  according  to  their  key 
characteristics,  and  to  identify  corrective  actions  permissible  within  the  project’s  existing 
authorization. 

The  DMS  has  three  components — the  Manual,  the  Data  Manager  (Johnston  and  Craig  2002),  and 
the  Analytical  Toolbox.  The  Analytical  Toolbox  is  a  collection  of  calculation  tools  and  guidance 
documents  assembled  for  understanding  and  quantifying  shoaling  conditions  and  for  formulating 
and  evaluating  alternatives  to  reduce  dredging.  The  calculation  tools  consist  primarily  of 
analytical  and  simple  numerical  models  for  estimating  hydrodynamics  and  sediment  transport  in 
and  around  channels  for  a  variety  of  situations.  These  tools  are  easy-to-use,  stand-alone 
applications. 

The  EbbJet  Calculator  estimates  tidal  current  velocity  along  an  entrance  channel  of  an  inlet  based 
on  the  theory  of  Ozsoy  and  Uni uata  (1982).  In  their  theory,  the  velocity  distribution  of  an  ebb- 
or  flood-tidal  jet  is  calculated  including  lateral  mixing  and  entrainment,  bottom  friction,  one¬ 
dimensional  bathymetric  changes,  and  ambient  currents.  The  EbbJet  Calculator  also  provides 
estimates  of  sediment  transport  through  equations  developed  by  van  Rijn  (1984a,  1984b,  and 
1984c).  Knowledge  of  the  flow  patterns  and  sediment  transport  around  inlets  will  aid  in 
optimizing  channel  alignment  and  estimating  shoaling  rates  under  possible  planned 
modifications. 
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THEORY:  The  EbbJet  program  employs  the  theory  of  Ozsoy  and  Unliiata  (1982)  to  simulate 
tidal  inlet  hydraulics  on  ebb  tide.1  Analytic  equations  are  developed  from  the  depth-averaged 
equations  of  motion  taking  into  account  lateral  mixing  and  entrainment,  bottom  friction,  and  one¬ 
dimensional  bathymetry  change.  The  EbbJet  program  does  not  take  into  account  the  ambient 
current.  The  theory  is  based  on  the  assumption  of  an  idealized  inlet  with  equal  shoreline  offsets. 
The  offshore  bathymetry  consists  of  straight  and  parallel  contours.  The  Darcy-Weisbach  friction 
factor  is  used  to  model  bottom  shear  stress  and  remains  constant  over  the  solution  domain. 

Consider  the  turbulent  jet  flow  produced  by  water  ebbing  from  an  inlet  into  a  semiinfinite  ocean 
as  depicted  in  Figure  1.  In  the  figure,  x  is  the  offshore  distance,  y  is  the  longshore  distance  from 
the  center  of  the  inlet,  b(x)  is  the  half-width  of  the  jet,  r(x)  is  the  half- width  of  the  jet  core 
(region  where  the  velocity  is  the  same  as  the  inlet  entrance),  w(x, y)  is  the  x  component  of 
velocity  at  any  point  in  the  jet,  wc(x)  is  the  x-component  of  velocity  along  the  jet  center  line,  and 
h{x)  is  the  depth  (from  still-water  conditions).  The  constants  ho,  bo,  and  u0  describe  the  depth  at 
the  inlet,  the  half-width  of  the  inlet,  and  the  velocity  at  the  inlet. 

The  reduced  equations  of  motion  (conservation  of  mass  and  momentum)  take  the  following 
form: 


dhu 1  dhuv  f  o  Id,, 

dx  dy  8  p  dy  y 

dhu  dhv 

- + - =  0 

dx  dy 


(1) 


where  h  is  the  local  water  depth,  u  and  v  are  the  x-  and  y-  directed  velocities,  respectively,/^  the 
Darcy-Weisbach  friction  factor,  p  is  the  density  of  water,  and  Fxy  is  the  cross  coupling  tenn  of 
the  depth-integrated  stress  tensor  (including  turbulent  fluctuations). 

To  solve  these  equations,  the  velocity  distribution  is  assumed  to  be  self-similar  with  respect  to 
the  nondimensional  coordinate: 


<r= 


6(x) 


(2) 


The  similarity  profile  is  adopted  from  Abramovich  (1963).  It  has  the  fonn 

‘jpp=J& 

*k(x) 


(3) 


1  The  EbbJet  Calculator  is  also  capable  of  estimating  the  current  velocity  and  potential  sediment  transport  for  the 
flood  current.  For  the  purposes  of  this  CHETN,  only  calculation  of  the  ebb  current  is  described. 
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where 


m 


0  ;  1<§ 
(1-f5);  0<^<1 
i  ;  \<o 


l=kli± 

I-/-  /  b 


(4) 


Zone  of  Flow  Establishment  (ZOFE)  —  0<  x  <  xs 
Zone  of  Established  Flow  (ZOEF)  —  x  >  xs 

Two  regions  within  the  jet  cross  section  characterize  the  ZOFE:  a  core  region  where  the  velocity 
is  the  same  as  the  velocity  issuing  from  the  orifice  (inlet,  x  =  0)  and  a  turbulent  mixing  region 
where  the  velocity  is  less  than  that  along  the  jet  center  line.  The  turbulent  mixing  zone,  having 
zero  width  at  the  jet  orifice,  grows  to  cover  the  entire  jet  at  the  end  of  the  ZOFE.  Conversely,  the 
core  region  extends  through  the  entire  jet  cross  section  at  the  orifice  and  reduces  to  zero  at  the 
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end  of  the  ZOFE.  The  ZOEF  is  characterized  by  a  turbulent  mixing  region  that  extends 
throughout  the  jet  cross  section.  In  the  ZOEF,  the  jet  is  fully  developed  and  the  core  region  has  a 
width  of  zero. 


Characteristics  of  each  of  these  regions  allow  for  simplification  of  the  differential  equations  to  an 
analytical  solution  given  the  following  boundary  conditions: 


m=b0 

r(0)  =  b0 
uc{0)  =  u0 
r(x)  =  0  ;  x>x5 

Uq  ^  X  ^  —  Wq  ,  X  <  X^ 

v[x,6(x)]  =  auc  (x) 


(5) 


The  final  condition  states  that  the  longshore  velocity  (entrainment  velocity)  at  the  jet’s  edge  is 
proportional  to  the  velocity  along  the  jet  center  line  through  the  constant,  a,  where  a  is  equal  to 
0.036  in  the  ZOFE  and  0.050  in  the  ZOEF  (empirically  determined). 

Defining  the  following  dimensionless  parameter  allows  for  reduction  of  the  equations  of  motion: 


S  =  b  =  =  R(Q  =  {-,  m  =  ^~,  = 

b0  8  fiQ  h0  b0  b0  u 


(6) 


The  equations  of  motion  are  then  solved  to  yield  the  following  results  for  the  ZOFE: 


TT  ,  D  0.45F-0.314G  0.683G-0.55F 

U  —  1 ,  K  = - ,  B  —  - 


0.134  H 


0.134  H 


(7) 


where 


F  =  exp 


It fL  ■ 


0  b  J 


In  the  ZOEF,  the  following  equations  apply: 


F  L 

R  =  0,  U  =  -=,  B  = - 

VZ  0.316F/F 


(8) 


(9) 
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where 


L=^r\ ^  =t~  (io) 

0.45  J  b0 

Both  sets  of  equations  yield  solutions  for  b,  r,  and  uc  (and  hence  w(x,y))  given  the  inputs  of  bo, 
wo,/,  and  h(x). 

Sediment  transport  is  estimated  from  the  equations  developed  by  van  Rijn  (1984a,  1984b,  and 
1984c)  and  the  current  velocities  calculated  within  the  EbbJet  Calculator.  The  total  sediment 
transport  rate  includes  bed  load  and  the  suspended  load.  Van  Rijn  computed  bed-load  sediment 
transport  as  the  product  of  the  saltation  height,  the  particle  velocity,  and  the  bed-load 
concentration.  Saltation  height  and  particle  velocity,  for  a  variety  of  sediments  and  flow 
parameters,  is  solved  numerically  with  the  equations  of  motion.  The  model  calibrated  with 
laboratory  data.  Measurements  of  bed-load  transport  provided  empirical  determination  of  the 
sediment  concentrations  in  the  boundary  layer. 

Before  calculating  the  bed-load  transport,  the  dimensionless  particle  size  parameter  and  the 
transport  stage  parameter  need  to  be  defined.  The  particle  size  parameter,  D*  is  found  using  the 
following  equation: 


I)<  —  /) 


'50 


(s-l)g 


1/3 


(ID 


In  the  equation,  D50  is  the  median  diameter  of  the  bed  material,  5  is  the  specific  gravity  of  the 
sediment,  g  the  acceleration  of  gravity,  and  v  is  the  kinematic  viscosity  of  the  fluid.  The 
nondimensional  transport  stage  parameter,  T,  is  calculated  via  the  following  equation: 

•\2_( 

I  I  U%,  Qy 

( u*,cr ) 


(12) 


where  u'*  is  the  bed  shear  velocity  related  to  grain  size  (on  the  bed)  and  u*cr  is  the  critical  shear 
velocity  according  to  Shields'  Curve.  The  parameter  u'*  is  an  estimate  for  the  average  shear 
velocity  at  the  upsloping  part  of  the  bed  forms. 


Given  D*  and  T,  calculation  of  the  bed-load  transport,  qt,  is: 

-7,2.1  ^ 


%=»502 


0.053 


A03 


V(s_1)k 


(13) 


Bed-load  transport  is  given  in  units  of  volume  of  transported  material  per  unit  time  per  unit  width 
(van  Rijn  1984b). 
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Van  Rijn  (1984b)  computed  the  depth-integrated  suspended  sediment  transport  rate  as  the 
product  of  local  concentration  and  flow  velocity.  A  reference  concentration  is  derived  from  the 
equations  for  bed-load  transport.  Measured  concentrations  provided  the  calibration  for  the 
equations. 

The  following  equations  result  from  the  integration  of  the  differential  equation  describing  the 
vertical  distribution  of  sediment  concentration  in  a  steady  uniform  current: 


<ls 


_  u*cg 


K 


a 

h—a 


Z' 


0.5  h 


h  —  z ^ 


Z' 


In 


Vzo) 


dz+  j  e-4Z'(^-0.5)ln 


0.5  h 


f  \ 

z 

dz 

yZ0) 

(14) 


where  u*  is  the  friction  velocity,  ca  is  the  reference  concentration,  a.' is  von  Kannan’s  constant,  a 
is  the  reference  concentration  elevation,  h  is  the  local  water  depth,  Z’  is  a  suspension  parameter, 
z o  is  the  zero  velocity  level  above  the  bed,  and  z  is  the  vertical  coordinate  (positive  upwards  from 
the  bed).  This  equation  is  represented  by  the  simpler  expression: 


q  =  Fsuhca 


Fs=- 


(  A2' 
'  a  ' 


ynj 


(  A1-2 
'  a  ' 


X 


1-- 

h 


(1.2-  Z') 


(15) 


for  the  range  0.03  <  Z  <  3  and  0.01  <  a/h  <0.1  with  a  maximum  error  of  approximately 
25  percent.  In  the  equations,  the  reference  concentration  is  defined  as 


0.035  T>5q  Tls 
2.3  a  £>0-3 


(16) 


To  avoid  errors  associated  with  setting  the  reference  elevation,  a,  too  low,  van  Rijn  (1984b) 
assumes  a  reference  elevation  equal  to  one -half  the  bed  fonn  height  or  the  equivalent  roughness 
height,  if  the  bed  form  dimensions  are  unknown.  The  minimum  reference  elevation  should  be  no 
less  than  1/1 00th  the  local  water  depth. 

Finally,  the  suspension  parameter  in  the  suspended  load  equation,  Z’,  is  defined  by 

Z'  =  Z  +  (p  (17) 

where 


Z  = 


Pku* 


(18) 
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<p  =  2.5 


,0.8 


C0c 


f 

c_a_ 

yC0 


,0.4 


0.01<  — <1 

U* 


(19) 


P  =  1  +  2 


4,,  A2 

y  U*  / 


0.1<  — <1 

U* 


(20) 


In  the  previous  equations,  cqv  is  the  sediment  fall  velocity  and  Co  is  the  maximum  volumetric  bed 
concentration  (=0.65). 


EBBJET  CALCULATOR  INTERFACE:  The  EbbJet  Calculator  is  a  front-  and  back-end 
model  with  integrated  interface  developed  for  solving  the  equations  as  previously  described. 
Upon  starting  the  program,  the  user  is  presented  with  the  console  shown  in  Figure  2.  The 
console  consists  of  three  areas:  the  menu  bar,  the  graph  display  area,  and  the  graph  display 
option  buttons.  The  menu  bar  runs  the  calculations  for  the  program;  the  graph  display  area 
shows  the  results  of  the  calculations,  as  either  line  or  contour  graphs;  and  the  graph  display 
option  buttons  select  which  data  are  presented  and  how  they  are  presented. 
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Figure  2.  EbbJet  console 


To  begin  a  simulation,  the  user  can  either  load  a  simulation  from  an  external  input  file  or  enter 
the  information  manually.  To  enter  the  information,  select  the  menu  item  Calculate\lnlet 
Hydraulics.  This  action  causes  the  Inlet  Hydraulics  dialog  box  to  appear  (Figure  3).  Here  the 
user  can  either  load  information  from  a  previously  saved,  formatted  input  file  or  manually  enter 
the  data  requested.  Required  data  includes  the  depth  at  the  inlet,  the  inlet  half-width,  the  depth- 
averaged  velocity  at  the  inlet,  the  Darcy-Weisbach  friction  factor,  a  description  of  the  offshore 
area  (as  either  a  constant  slope  or  an  input  profile),  and  the  units  of  all  calculations  (metric  or 
English). 
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Figure  3.  Inlet  Hydraulics  dialog  box 

After  the  required  information  has  been  entered,  the  Calculate  Inlet  Hydraulics  button  will 
solve  the  hydraulics  equations.  The  program  returns  to  the  main  console  where  the  user  can 
view  the  results  either  as  line  graphs  (Figure  2)  or  as  contour  graphs  (Figure  4).  Calculation 
results  available  for  display  as  line  graphs  are  the  input  bottom  elevation,  the  jet  half-width,  core 
half-width,  center  line  velocity,  and  center  line  sediment  transport  (enabled  after  performing  the 
sediment  transport  calculations).  Contour  graphs  available  for  display  are  cross-shore  velocity 
magnitude  and  sediment  transport  magnitude  (enabled  after  performing  the  sediment  transport 
calculations). 


Figure  4.  Contour  graph  of  cross-shore  velocity  magnitude 
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The  program  is  structured  so  that  a  hydraulic  simulation  must  be  completed  before  the  sediment 
transport  calculations  can  proceed.  To  perform  sediment  transport  calculations,  the  user  selects 
the  menu  option  Calculate\lnlet  Sediment  Transport.  This  action  causes  the  Inlet  Sediment 
Transport  dialog  box  to  appear  (Figure  5).  Again,  the  user  can  either  load  the  requested 
information  from  a  previously  saved,  formatted  file  or  enter  the  information  manually.  Required 
sediment  information  includes  D50,  D90,  sediment  specific  gravity,  and  the  sediment 
distribution’s  geometric  standard  deviation.  Additionally,  the  fluid  density  and  kinematic 
viscosity  are  required  for  the  calculations.  The  fluid  properties  may  be  input  manually  or 
calculated  from  the  fluid  temperature  and  salinity. 


Figure  5.  Inlet  Sediment  Transport  dialog  box 

Upon  entering  the  requested  information,  the  user  presses  the  Calculate  Sediment  Transport 
button  to  perform  the  calculations.  The  program  then  returns  to  the  main  console  where  the 
results  may  be  displayed  as  either  line  graphs  or  contour  graphs  (Figure  6). 
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Figure  6.  Contour  graph  of  sediment  transport  contours 

In  addition  to  the  features  previously  described,  the  EbbJet  Calculator  contains  several  other 
capabilities  accessed  through  the  following  menu  items: 

a.  File\Save  Simulation  —  This  action  launches  a  Save  File  dialog  box.  The  saved  fde 
contains  only  the  inputs  entered  into  the  Inlet  Hydraulics  dialog  box  and  the  Inlet 
Sediment  Transport  dialog  box.  This  file  can  be  loaded  at  a  later  date  to  rerun 
simulations  and  view  results. 

b.  Edit\Copy  Chart  to  Clipboard  —  Selecting  this  option  copies  the  visible  graph  (line  or 
contour)  to  the  clipboard.  The  graph  is  then  available  for  pasting  into  any  document  or 
application  as  a  Windows  metafile. 

c.  File\Export  to  File  —  Selecting  this  menu  item  exports  the  temporary  files  used  to  create 
the  graphs  to  ASCII  files.  This  makes  the  results  of  the  calculations  available  for  display 
in  an  alternate  software  package. 

d.  File\Print  —  This  menu  item  displays  the  Print  dialog  box.  Selecting  the  appropriate 
printer  and  clicking  OK  will  print  the  graph  showing  in  the  Graph  Display  Area. 

e.  View\Variable  Summary  —  This  menu  option  displays  the  Variable  Summary  window 
which  contains  all  the  information  listed  in  the  Inlet  Hydraulics  dialog  box  and  the  Inlet 
Sediment  Transport  dialog  box.  Text  in  the  window  may  be  selected  and  copied  to  the 
clipboard  for  subsequent  pasting  into  a  text  editor. 

f.  Help  —  An  extensive  help  file  is  also  available  through  the  menu  bar. 

Calculating  the  currents  and  sediment  transport  at  a  tidal  entrance  with  the  EbbJet  Calculator 
provides  a  quick  overview  of  the  flow  characteristics  at  the  entrance.  This  information  can  help 
evaluate  channels  near  inlets.  For  example,  a  hydrographic  survey  can  be  designed  to  adequately 
cover  potential  shoaling  areas  due  to  tidal  sediment  transport.  By  entering  spring  tide  conditions 
into  the  EbbJet  Calculator,  areas  of  active  sediment  transport  can  be  identified.  The  results  will 
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provide  a  rough  estimate  of  the  offshore  extent  of  the  ebb  tidal  shoal.  If  the  results  resemble  the 
output  in  Figure  6,  the  project  manager  can  specify  that  the  survey  should  extend  approximately 
1,524  m  (5,000  ft)  seaward  of  the  inlet  entrance. 

CONCLUSIONS:  The  EbbJet  Calculator  allows  quick  estimation  of  the  ebb  or  flood  jet  issuing 
from  an  inlet  based  on  the  analytical  solution  of  flow  from  Ozsoy  and  Unluata  (1982).  Estimates 
of  sediment  transport  are  developed  from  the  sediment  transport  equations  published  by  van  Rijn 
(1984a,  1984b,  and  1984c).  Developed  for  inclusion  in  the  Diagnostic  Modeling  System 
Analytical  Toolbox,  the  EbbJet  Calculator  provides  estimates  of  flow  velocity  and  associated 
sediment  transport  at  inlets  that  will  aid  in  optimizing  channel  location  and  estimation  of 
shoaling  rates. 

OBTAINING  SOFTWARE:  The  EbbJet  Calculator  is  distributed  from  the  DMS  Web  site 
located  at  http://www.taylorengineering.com/DMShome/DMSDefault.htm  and  at  the  Coastal 
Inlets  Research  Program  Web  site  located  at  http://cirp.wes.army.mil/cirp/cirp.html,  then  go  to 
Products  and  Tools. 

ADDITIONAL  INFORMATION:  Questions  about  this  CFIETN  may  be  addressed  to 
Dr.  Mark  S.  Gosselin  (904-731-7040;  e-mail:  mark@oea-inc.com ).  The  authors  appreciate 
review  of  this  CHETN  by  colleague  Dr.  R.  Bruce  Taylor  and  by  Dr.  Nicholas  C.  Kraus,  Principal 
Investigator  of  the  Diagnostic  Modeling  Work  Unit,  U.S.  Army  Engineer  Research  and 
Development  Center,  Coastal  and  Hydraulics  Laboratory  (601-634-2016;  e-mail: 
Nicholas.  C. Kraus@erdc.usace.army.mil).  This  CHETN  should  be  cited  as  follows: 

Gosselin,  M.  S.,  and  Craig,  K.  R.  (2002).  “Estimating  along-channel  flow  and 
sediment  transport  at  tidal  entrances  —  EbbJet  Calculator,”  Coastal  and  Hydraulics 
Engineering  Technical  Note  CHETN-IV-48,  U.S.  Army  Engineer  Research  and 
Development  Center,  Vicksburg,  MS. 

http://chl.wes.army.mil/library/publications/chetn/ 
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